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We explore in detail a dilatonic scalar-tensor theory of gravity inspired by large extra dimensions, 
where a radion field from compact extra dimensions gives rise to quintessence in our 4-dimensional 
world. We show that the model can give rise to other types of cosmologies as well, some more 
akin to fc-essence and possibly variants of phantom dark energy. In our model the field (or radius) 
stabilization arises from quantum corrections to the effective 4D Ricci scalar. We then show that 
various constraints nearly determine the model parameters, and give an example of a quintessence- 
type cosmology consistent with observations. We show that the upcoming SNAP-experiment would 
easily distinguish the present model from a constant A model with an equal amount of dark energy, 
but that the SNAP-data alone will not be able distinguish it from a A model with about 5% less 
dark energy. 

PACS numbers: 98.80.-k, 98.80.Cq, 95.35.+d, 04.50.+h 



I. INTRODUCTION 

The discovery that the expansion of the universe is ac- 
celerating implies the existence of an extremely small, 
but non-zero vacuum energy density. In summary, the 
observations have led to a new cosmological concor- 
dance model, according to which fi tot = 1.02 ± 0.02, 
57a ~ 0.73, and Q m w 0.27. That is, the universe is 
flat and dominated by a "dark energy" component Q\ 
which has an equation of state close to that of a bare 
cosmological constant today. The physical nature of the 
dark energy component is not known, but according to 
the popular quintessence scenario 0, 0, H , it corresponds 
to energy stored into a classical background scalar field 
with some very particular properties. An interesting class 
of quintessence models can be characterized by a La- 
grangian [IS HQ 
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where V p (x) is some power-law modification to the dom- 
inant exponential factor. A desirable feature of such ex- 
ponential potentials is that a small change in the field 
value x can cause a large change in the energy density 
p x . Moreover, they often allow scaling solutions, where 
the kinetic and potential energies of x maintain a fixed ra- 
tio, leading to a constant equation-of-state parameter w x 
and to p x scaling exactly as a power of the scale factor. 
This in particular makes it possible for the quintessence 
field to mimic a cosmological constant. 

Generic quintessence scenarios often rely on an ad. hoc. 
quintessence field with some suitably adjusted poten- 
tial. A more complete physical theory should explain the 
origin of the — possibly effective — quintessence field 



and naturally predict the exponential form of its poten- 
tial. Remarkably, these requirements can be met in a 
model with large extra dimensions, studied by Albrecht, 
Burgess, Ravndall, and Skordis (ABRS) 0,11 El- In the 
ABRS scenario the quintessence field x corresponds to 
the radion mode, which after a compactification of the 
originally 6-dimensional space-time manifold shows up 
as an effective Jordan-Brans-Dicke field in the residual 
4-dimensional action. A crucial ingredient of the ABRS 
scenario is its introduction of a scalar field <f> living in 
the 6D bulk space, causing two key features that make 
the model special: First, the zero-point fluctuation of the 
bulk scalar field <j> induce a Casimir potential Uq ~ Uq/ r 4 
for the radion, which eventually acts as an effective cos- 
mological constant. Second, the bulk scalar field 4> also 
induce logarithmic corrections to both kinetic terms and 
the radion potential 0, EH, which are crucial for the 
radius stabilization mechanism. The proper amount of 
dark energy can arise in a parametrically natural way 
from the Casimir potential Uq ■ In a model with d equally 
large extra dimensions, the effective strength of gravity is 



Hence, since the naturalness condition Uq ~ 1 
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requires a radius r ~ (10 _3 cV) _1 to yield the desired 
amount of dark energy, a 6D model corresponds to elec- 
troweak scale Mb ~ 1 TeV, whereas for d = 3 one would 
need a lower scale Mb ~ 1 GeV. Seen the other way 
around, setting the fundamental scale to TeV in a model 
with d > 2, would require fine-tuning Uq to extremely 
small values. In this sense the 6D case is indeed special, 
since it naturally yields both the electroweak scale and 
the desired amount of dark energy. After obtaining the 
general form of the action in section [HI we will set d = 2 
for the remainder of the paper. 

We will in this paper explore a class of dilatonic scalar- 
tensor theories of gravity (STG models for short) inspired 
by the ABRS scenario. The above extra-dimensional 
models can indeed be recast as a dilatonic STG model 
by simply defining the field ip = logMbr. Similar to the- 
ories with large extra dimensions, a dilatonic factor e Av 
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allows the fundamental scale of gravity Mb to be much 
lower than Mp\. Moreover, when written in the Einstein 
frame, any dilatonic STG model will introduce an ex- 
ponential behaviour to the potential of the theory. In 
a general scalar-tensor theory, not only the scalar field 
potential can cause an observed acceleration of the uni- 
verse, but also the functions associated with the kinetic 
terms of the theory. This paper will however focus on 
the possibility of an accelerated expansion originating 
in the potential term. Although the connection to the 
ABRS scenario provides additional motivation, it should 
be stressed that our results are independent of origin and 
the presented scalar-tensor theory can indeed be viewed 
as the starting point. 

The present model is strongly confronted by the exist- 
ing astrophysical constraints, and most notably by preci- 
sion tests of General Relativity in the solar system [l2| . 
GR-constraints were already considered by the authors 
of ref. [8|, but our results differ from theirs; in fact we 
found that they rule out the original scenario when cor- 
rectly accounted for. The problem arises from the fact 
that the perturbative loop corrections coming from the 
bulk scalar field </>, which were treated as small quanti- 
ties in ref. , actually must be of order unity to fulfil the 
observational constraints. This is the case in particular 
for the scaling function A(ip) for the effective 4D Ricci 
scalar R g in the Jordan frame: At tree level A(ip) = 1 
whereas the GR-constraint require that A((p) sa today. 
Although this poses no problem when one takes the STG 
model of Eqn. J2J) as a pure phenomenological starting 
point, we point out that it is indeed also possible to fulfil 
this requirement in the context of a converging pertur- 
bation theory within the extra dimensional motivation. 
This is technically possible if the loop corrections come 
from many independent very weakly coupled bulk fields 
<pi . The sacrifice in this scheme is that the original radius 
stabilization mechanism based on the Casimir potential 
does not work. 

However, we have found that our model allows a new 
stabilization mechanism, based on the (^-scaling of the ef- 
fective 4D Ricci scalar. Indeed, suppose that the function 
A(ip) has a zero or becomes very small for some ip — ip*. 
It turns out that a root A((p*) = shows up as a sin- 
gularity in the effective potential in the Einstein frame, 
that confines the field (and hence the radius in the extra- 
dimensional context) to the range <p < tp*. We shall see 
that this mechanism works in the simplest perturbative 
parametrization of the loop-corrected action. However, 
the mechansim is much more generic; we can state that 
in any model where A(ip) has a zero, for some large value 
ip — if*, the effective potential becomes confining and de- 
velops a minimum close to ip„, where A(tp m i n ) « 0. If 
strange at first sight, it should be noted that this scal- 
ing arises naturally in the present model, and it is in 
fact all but a necessary condition for satisfying the GR- 
constraint at present. Conversely, in the context of this 
stabilation scheme the GR-constraint becomes almost au- 
tomatically satisfied. 



Astrophysical and theoretical constraints nearly fix the 
relevant model parameters. Moreover, for a given set of 
parameters, there is practically a unique set of initial 
conditions that can give rise to correct values for the 
cosmological parameters today. We present a particu- 
lar quintessence-like solution which obeys all constraints, 
starting from time before nucleosynthesis until today. We 
then argue that our model makes a generic prediction for 
the evolution of the equation-of-state parameter w to t (z) , 
which in principle makes it distinguishable from a flat, 
constant A cosmology. We then show that the sensitivity 
of the forthcoming SNAP-experiment 01 should be good 
enough to separate our model from a constant A model 
with equal amount of dark energy today. However, the 
SNAP-data 0] alone will not be able to break a degen- 
eracy between our model and a A model with slightly 
different amount of dark energy (about 5% less). 

While our scalar-tensor action is very simple in the 
Jordan frame, it introduces a very complicated structure 
for the kinetic term of the scalar field ip when written in 
the Einstein frame. Even with the little room one has to 
play after applying the constraints, it is in fact possible to 
arrange model parameters such that this term becomes 
very large and possibly singular, or even negative. In ei- 
ther of the former cases the model classifies rather as a 
fc-essence than a quintessence scenario, and in the latter 
it gives rise to variants of recently introduced phantom 
cosmologies [l5j . In this paper we mostly concentrate on 
non-phantom models. While we do give some examples 
of parametrizations leading to phantom cosmologies, it is 
not yet clear to us whether it is possible to make a transi- 
tion into such a phase and make the scenario compatible 
with observations. 

The paper is organized as follows. We begin in sec- 
tion [H] by introducing the scalar-tensor theory, and show 
its motivation within extra-dimensional theories, by re- 
viewing the basic ABRS scenario and its reduction to an 
effective 4D gravity theory. We continue in section lHll bv 
discussing how the effective 4D theory relates to other 
dark energy models, including fc-essence, quintessence, 
and phantom cosmology. We also derive an effective ac- 
tion for the model in terms of the field x m t he Einstein 
frame, and from it the coupled cosmological evolution 
equations for \ an d various density components. Section 
HVI beeins with a discussion of a problem with the pertur- 
bative expansion employed by ABRS and then outlines 
a perturbatively reliable scheme with a new mechanism 
for the radius stabilization. Having obtained a consistent 
scenario, we continue in section^by deriving constraints 
on the model parameters, arising from various theoret- 
ical and astrophysical considerations. In section IVII we 
then present a realistic cosmological evolution from times 
before nucleosynthesis until present that satisfies all ob- 
servational constraints. We also compute the prediction 
of the model for the supernovae magnitude-redshift rela- 
tion and compare it with cosmological constant models. 
Finally, section IVHI contains our summary and the con- 
clusions. 
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II. THE MODEL 



described by an effective action 



The aim of this paper is to explore the alternative cos- 
mologies of a class of scalar-tensor theories of gravity 



J 



S 



d 4 x, 



-*^-(A{y)R s -2B{<p){dyy 



- U(<p) + C n 



(2) 



r 



Here Mb is the fundamental scale of gravity, R g is the 
Ricci scalar of the metric g^, and £ m is the classsical 
Lagrangian of all matter contained within the universe. 
The dilatonic factor e x,p allows the scale Mb to be much 
lower than Mpi, but apart from this behaviour, A(ip) 
and B(tp) can be arbitrary functions. We will however 
see that the theory gives rise to interesting cosmologies 
already for the simple case when A(ip) and B(ip) are lin- 
ear in ip. One of the two functions A(tp) and B(<p) can of 
course be set to 1 by a redefinition of tp, but we have kept 
the above form since it will make the connection to extra- 
dimensional scenarios more explicit. While the detailed 
form of U(ip) is not constrained in a generic STG model, 
we will set U(ip) oc e _4v , since this will connect our class 
of scalar-tensor theories to the brane-world scenario to 
be discussed below. Note that for A(ip) = B(<p) = 1 and 
U(<p) — the action corresponds to a basic Jordan- 
Brans-Dicke theory with coupling parameter u> = 2/ A 2 . 

The above scalar-tensor theory, Eqn. (J2J, can be 
viewed as a starting point of this paper, but is should 
be stressed that the form of the action is motivated by 
extra-dimensional scenarios 0, @ ■ Because this will also 
lead to a more specific prediction of the form of the func- 
tions A, B and U, the present model is potentially more 
predictive than previous studies of scalar-tensor gravity 
and quintessence, where the scalar field corresponded to 
a dilaton from string theory 0, 0, or the origin of the 
field was not considered in depth |l8l ITsj. Also, in the 
present model, all mass scales and parameters have per- 
fectly natural values and no fine-tuning is needed for end- 
ing up at the desired amount of dark energy at present. 



A. The Brane World 

Let us consider a brane world scenario in which our 
(3 + l)-dimensional space is a brane embedded in a 
(4 + d)-dimensional bulk space. All standard model fields 
are confined to the four-dimensional brane while gravity 
propagates in all dimensions. It is well known that if 
the size of the extra dimensions r is large, the scale of 
physics on the brane Mb can be of electroweak scale, 
hence provi ding a compelling solution to the hierarchy 
problem [2(3, |2l| . The effective 4D gravitational coupling 
in this scenario is (87rG c ff ) _1/2 = M h (M h r) d / 2 ~ M P i ~ 



10 18 GeV, corresponding to a very large dimensionless 
parameter M b r - 10 30 / d TeV/M b . 

The full action of the theory can be divided into two 
terms, S = Sb + 5b, corresponding to the bulk and the 
brane contributions respectively. The leading part of the 
bulk action is 



S B 



M 



2+d 



(3) 



where 1Z is the Ricci scalar of the (4+<i)-dimensional met- 
ric Gmn and the ellipses represent possible bulk fields. 
The detailed form of the brane action is not important; 
the crucial part is that it only depends on the four- 
dimensional metric g^ v : 



d 4 s 



gC m (g^, ...) 



(4) 



All matter fields, denoted by ellipses in £ m , are trapped 
on the brane and do not carry stress-energy into the extra 
dimensions. The bare cosmological constant A is here 
assumed to vanish; at present our model does not offer 
any insight as to how this is accomplished. Further, we 
assume that the (4 + e?)-dimensional metric has a block- 
diagonal structure: 



JMN 



9»v{%) 
Q 2 {x)h mn (y) 



(5) 



That is, the (4 + rf)-dimcnsional space-time manifold con- 
sists of a four-dimensional brane with metric g^ Ul con- 
nected to a d-dimensional surface with metric h mn , such 
that the size of the extra dimensions, g(x) = M\ 1 r(x), is 
allowed to vary as a function of the position on the brane. 



B. Effective 4D Theory 

We are obviously interested in the effective four di- 
mensional theory on the brane at large scales. It can 
be obtained from Eqns. 100 via standard dimensional 
reduction. Indeed, integrating out the extra spatial di- 
mensions at the tree-level one finds the effective action: 
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SeS = 



M 2 (M b r) d 



Rg - d(d-l) 



Or 



{M h rf 



d d y^{-) d hR h + C m (g) 



(6) 



where Rg and Rh represent the Ricci scalars of g^ v and 
hmn , respectively, and the normalization of h mn is chosen 
such that M d J d d yy / (-) d h = 1. 

The action © is a modification of the regular Ein- 
stein Hilbert theory of gravity, and it can be recast into 
a dilatonic scalar-tensor model by a change of variables 
if = log M b r. It is well known that the presence of such 
a scalar field will yield long range forces in addition to 
the regular gravitational force originating in g^ v . As we 
will see in section precision measurements of general 
relativity on the solar system scale impose a strong con- 
straint a ps 0, where a is the coupling strength of the 
scalar field to matter (see Eqn. (|28[) below). In the tree 
level action 10 we would obtain a 2 — d/(d + 2), so that 
a ~ 1 at tree level in a generic model with d compactified 
large extra dimensions. 

It is therefore essential that the action can be modi- 
fied from its tree level form. If for example (M b r) d — > 
(M h r) d A, then a can be small if dA/dtp + dA sa for 
large M b r, where ip = log M b r. A simple example of 
this is a case when A is a linear function with a zero 
in log M b r. Remarkably, such logarithms are generically 



induced by quantum corrections from bulk fields fa liv- 
ing in large extra dimensions 0, El El ■ These include 
corrections to the kinetic terms, in particular an A-term 
as requested, but also a Casimir potential Uc, which for 
toroidal compactification has the leading term 



is, -5 



The Casimir potential also receives logaritmic quantum 
corrections. 

The curvature term proportional to Rh in Eqn. © 
is non- vanishing except for Ricci flat compactifications. 
Indeed, if Rh ~ ^ ms t erm would be too large 

to be consistent with the observed amount of the dark 
energy and the value of Newtons constant for any number 
of extra dimensions. We are therefore led to consider 
only topologies where the contribution from Rh vanish, 
such as the trivial case when the extra dimensions are 
Ricci flat (Rh = 0) (ijj. These considerations lead us 
to parametrize the action of the effective 4D theory as 
follows: 



J 



SeS = 



d 4 x 



M b 2 (M b r)° 



A{r)R g - B{r)d{d 




r 



(8) 



where 



A(r) 
B(r) 
C{r) 



1 + a log M h r . 
l + 61ogM b r, 
1 + c log Mbr , 



(9) 
(10) 
(11) 



and a, 6, and c are unknown but small parameters pro- 
portional to some dimensionless coupling constants of the 
theory 

Let us now return to the original scalar-tensor theory. 
Indeed, comparing the action with the effective 4D 
action of the brane- world scenario, Eqn. (jSJ, we see that 
they are completely equivalent, and can be put into iden- 
tical form via the definitions 



if = log Mbr . 



1%) = C(f)U e-^ , 



(12) 



(13) 



and A(tp) = l + atp, B(tp) = 1 + bip, C(ip) — 1 + cip. Here, 
A = d and we now fix the number of extra dimensions 
to d = 2. For this particular choice of d, the naturalness 
condition Uq ~ 1 and the desire to obtain the observed 
amount of dark energy corresponds to the electroweak 
scale Mf, ~ 1 TeV. Indeed, this explicit form of the action 
will be the focus for the remainder of the paper. Note 
that we treat all fields in the action (JSJ), and hence ip, as 
purely classical fields. The corrections in the form of a 
Casimir potential and the coefficients A, B, and C arise 
from the presence of quantized bulk fields fa. 

The parametrization of the model JSJ was first derived 
by Albrecht, Burgess, Ravndal and Skordis, in the con- 
text of a six-dimensional model with a bulk scalar field 
<t> with renormizable cubic self-interaction and a toroidal 
compactification 0, • It has later been realized [24| that 
logarithmic corrections actually do not arise for this sim- 
ple toroidal compactification and boundary conditions. 
However, logarithmic corrections from bulk fields fa have 



5 



been found to be a generic outcome for a variety of com- 
pactification schemes |nj and moreover, it is likely 
that they are generated also for torii given twisted bound- 
ary conditions and/or orbifold symmetries |2^|- We thus 
believe that the effective actions and JSJ) have strong 
motivation based on large extra dimensional theories. 
Moreover, note that we basically derived the form © 
out of necessity: In a generic model with large extra di- 
mensions a large curvature term cx Rh must be absent 
and the constraints from general relativity force us to in- 
troduce a scaling function A such that A ss for large 
Mb?". The action © with the parametrizations (|51 lllf> is 
perhaps the simplest one compatible with these requests, 
but as will be further discussed in section ITTll the pre- 
cise forms of the scaling functions are not important; it 
is enough if A has a zero or just becomes very small for 
large M-^r, i.e. the present value of ip. 



C. Casimir Radius Stabilization 

In the context of large extra dimensions (LED), the 
original ABRS argument for large radius stabilization of 
the compactified dimensions was based on the particular 
form of the quantum corrected Casimir potential [EH EH 
in Eqn. ©: 



III. QUINTESSENCE, fc-ESSENCE OR 
PHANTOM ENERGY? 



We have so far written the model in the so called Jor- 
dan frame, where the metric g^ v is not canonically nor- 
malized. However, it is in the Jordan frame where all 
observations have their usual interpretation, since it is in 
this frame where the matter Lagrangian couples to the 
metric only and hence making it possible to define stan- 
dard clocks and rods that do not depend on the value 
of the scalar field. Nevertheless, the equations of motion 
are far more transparent in the Einstein frame (denoted 
below by a metric g^ v without a hat), where the gravi- 
tational action has the more familiar form: 



M?Ae Xv 



-Rg = 



Ml 

—£-Rq + additional terms , 



(16) 



where A = d and R g is the Ricci scalar of the Einstein 
metric g^ v . Let us define the commonly used variables 



F(ip) = Ae dip , Z(tp) = -d{d - l)Be 



dtp 



(17) 



where, F(ip) and Z((p) are the dimcnsionless functions 
multiplying the kinetic terms Rg and (dip) 2 , respectively, 
in the Jordan frame action (jHJ. Now, since schematically 
R g cx <7~J, Eqn. (|16|) implies the conformal transforma- 
tion 



U c (r) = U Q M£ 



1 



M h r 



[1 + c log M b r] 



(14) 



Uc{r) always has a runaway minimum at r — > oo. How- 
ever, if c is negative it develops another minimum at 
finite, but exponentially large scale (in comparision with 
the microscopic scale I = M^ 1 ): 



M h r 



(15) 



Fixing Mbr to the desired value ~ 10 15 indeed sets the 
scale of radiative corrections to a small and perturba- 
tively natural value |c| ~ 1/50. This simple argument 
neglects the coupling between the radius and the 4D met- 
ric in Eqn. (JSJ however. In reality one cannot make this 
assumption, and it will turn out that the other radia- 
tive corrections, and the A-teim in particular, will play 
a crucial role for the stabilization in later developments. 

Although the simplified line of reasoning above does 
not hold for the stabilization mechanism itself, it does 
set the right scale for the parameters a, b, and c needed 
by a successfull cosmology. In the LED context the small 
scale of these parameters is natural since they arise as 
radiative corrections. One could imagine this to also be 
the case for the original scalar-tensor theory, Eqn. J5J, 
even if one chooses to neglect its connection to the brane- 
world scenario, and we will henceforth adopt 1/50 as the 
natural scale for a, b, and c. 



9nu = F(^p)g^ v 



(18) 



This recasts the action |(SJ| with our choice d — 2, corre- 
sponding to the ABRS scenario, as 



Scff — 



d A x, 



^R g + 1 -P(ip)(dpf 



Ve(<p) + C m (g,ip) 



, (19) 



where 



fc 2 (y) 

Ml 



dlogFQg) 
dtp 



aip 



Z{<P) 
F{ip) 

2(1 + bip) 



1 



and 



F*{tp) 



= M^U 



1 



cip 



(l + aip) 2 



aip 



-Sip 



(20) 



(21) 



As seen from Eqn. I|21|l . any dilatonic theory with a non- 
zero U(ip) will give rise to an exponential potential in the 
Einstein frame. Note that we have not expanded any of 
the above expressions in perturbative parameters; it will 
later become clear why this indeed is the sensible thing to 
do. Written in the form (|19f) our effective theory has ob- 
vious resemblance with some popular toy models studied 
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earlier in the literature. First, if the radiative corrections 
were neglected (or they are neglibly small so that effec- 
tively a, b, c ~ 0) the potential reduces to a pure exponen- 
tial and the kinetic energy becomes canonical apart from 
a constant scaling. In this limit our theory thus becomes 
very similar to the original quintessence scenario |26j . 
We shall see later that the model, despite t does lead to 
quintessence-like cosmological solutions. 

In general, the function k 2 (p) gives rise to a non- 
canonically normalized kinetic term for <p, which is the 
trademark of another class of dark energy models called 
k-essence Hail]- In particular our model bears similar- 
ity to mixed fc-essence - quintessence models discussed in 
ref. [2]j. Observe in particular that in the Planck era, i.e. 
when tp ~ (corresponding to r ~ in the LED sce- 
nario) we have k 2 ((p) = 4M^ and V E (<p) = U M£. Thus 
the theory (|19|l is natural (contains no arbitrarily small 
parameters) if only Uq ~ 1. This is of course required in- 
dependently by the naturalness of our underlying model 
from a field theory point of view, and we shall see later 
that it does give rise to acceptable cosmologies as well. 

It is conceivable that for certain values of the parame- 
ters a and b, the function k 2 (tp) becomes negative. In this 
case our model corresponds neither to fc-essence nor to 
quintessence, but instead exhibits the features of phan- 
tom cosmologies |l5j . with an equation-of-state parame- 
ter w Ln < — 1 ■ (Of course, as was also pointed out in 
ref. (23, the phantom models could be viewed as form- 
ing a particular class of fc-essence models.) For exam- 
ple, one could imagine that our model initially starts as 
quintessence or fc-essence, and later in time changes its 
behaviour to a phantom cosmology, or vice versa. We 
shall show later such scenarios are indeed supported by 
natural values of the model parameters, but it is not yet 
clear to us whether they can be arranged dynamically 
in a model compatible with observations. For the major 
part of the paper we will however concentrate on the case 
where k 2 (p) > 0. 

Let us finally note that the conformal transformation 
(|18J) between the Einstein and Jordan frame metrics is 
singular at A{p st ) — 0, and that changing the sign of 
A will cause gravitons to carry negative energy in the 
Jordan frame and effectively reverse the signature of the 
Einstein metric. Perhaps fortunately this consideration 
appear mostly academic, as we shall see that it is practi- 
cally impossible for any classical evolution starting from 
a if smaller that tp* to move past the singularity, hence 
keeping the field confined to values below tp* at all times. 



term (in a pure fc-essence the potential is actually vanish- 
ing |2S|). whereas in quintessence the desired behaviour 
follows from the special form of the potential. We will 
find that our model can at times behave as quintessence 
and at others as fc-essence, as well as show some charac- 
teristics not typical to either of them. Indeed, introduc- 



ing 



dK(w) 

X = K{<p) with k(tp) = ' (22) 



dip 



we get the effective action 



S c g = 



M 2 



-R a 



;(dx) 2 



-V{ x )+£m(g,x) 



(23) 



where V( X ) = V E (p) = V B {K-\ X ))^ 

We have found an analytic expression for X m terms of 
ip from Eqns. (|20|) and l|22l) . It is rather complicated how- 
ever, and we defer the explicit solution to the appendix 
A. We plot the result for X as a function of ip in Fig. QJi). 
The dependence is smooth and linear except at a point 
where the derivative dxjd<p becomes very large (in fact 
infinite). As a result, the potential V is roughly expo- 
nential also as a function of x away from the singularity, 
i.e. automatically of the desired type for quintessence 
(Eqn. Q). We plot V( X ) for a representative choice of 
model parameters in Fig. Qd). V{x) is indeed exponen- 
tial for x —12. Moreover it has a minimum that leads 
to a stabilization of the field, which here corresponds to 
tp w 36 (or, in the LED context, M h r w 3.5 x 10 15 ). 
Note that the potential minimum is situated close to the 
dx/d^s-singularity in tp. This is so because the p-scsde is 
very strongly condensed with respect to x near the sin- 
gularity. We shall return to this and the other features 
seen in the potential function shortly. 



B. The Evolution Equations 

The action (|23|) . written in terms of Einstein frame 
variables, forms the basis of the cosmology to be studied 
in the rest of the paper. However, for a homogeneus 
and isotropic universe, it is the Jordan frame metric g^ v 
which is described by the Friedmann-Robertson- Walker 
line element: 



A. Non-Phantom Cosmologies 

When k 2 (tp) > we can make a further change of vari- 
ables such that the kinetic term has the standard canon- 
ical form [23 • In this case our theory belongs roughly to 
the class of fc-essence and quintessence models. The dif- 
ference between the latter two is that in A:-essence the be- 
haviour of the field is driven by its non-canonical kinetic 



dsj = g^dx^dx" = dt] - i?^j)dx 2 



(24) 



where ij and i?j are the measured Jordan frame time 
and scale factor, respectively. Since the relation between 
the Einstein and Jordan frame metrics is given by g^ v — 
F(p)g fll , 1 the universe is still of Friedmann-Robertson- 
Walker type in the Einstein frame, with time t and scale 
factor R given by: 



df 



Fdt, 



R = VFRj 



(25) 
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Fig. (a) 
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Fig. (b) 
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FIG. 1: (a) x ln terms of f>/ log 10 (= log 10 M^r) for the representative case a = —0.0276 and b = —0.031. (b) The potential 
V(x) for the same a and b with Uo = 2.60 and c = 0. The dashed line represents logF(x). Both x an d V(x) are given in 
Planck units. 



Hence, taking a variation with respect to x leads to an 
equation of motion |42j 

X + 3H± + V\x) + ^fJH = . ( 26 ) 

where H = R/R is the Einstein frame Hubble parame- 
ter. The explicit form of the derivative of the potential 
V'(x) = dV/dx is 



V\ X ) = M*U 



1 



cA - 2aC - 8 AC 



{M h rf 



A 2 JmA + a y -2AB 



, (27) 



and the coupling function a is given by 

dlogF- 1 ^ 
a(x) = V2M h 2- 

A + a/2 



1(2,4 + a) 2 -AB 



(28) 



These are all exact expressions (at the one-loop) and 
to evaluate them as a function of Xi one needs the 
inverse of the transformation 1)22(1 . i.e. for example 
A = 1 + aif _1 (x). Finally, T^" is the Einstein frame 
stress-energy tensor of radiation and matter, but since 
matter is pressureless and p = for radiation, we have 
simply T£ = p m . 

The Einstein frame Hubble parameter and the a-term 
couple the evolution of the field x to the evolution of 
the ordinary radiation and matter components. Varying 
the action l|23ll with respect to yields the Friedmann 
equations for the fiat, homogenous and isotropic cosmolo- 



gies: 



3M b 2 ff 2 = p. 



-2M^H 



P + P, 



(29) 
(30) 



where p and p are the total energy density and pressure 
in the Einstein frame: 



9 = Pn 
V = Pr 



F Pt + P x 
Px ■ 



(31) 
(32) 



At this point it is important to once again remind our- 
selves of the conformal transformation (JTSJ- In the Jor- 
dan frame, the strength of gravity is varying with the 
field <p and the model is characterized by an effective 
gravitational "constant" G c g: 



1 



M*F(<p) 



167rG cff 



(33) 



However, after conformal transformation to the Einstein 
frame, the apparent strength of gravity G is truly con- 
stant: 



1 



_ Ml 



16ttG 



(34) 



The dependence on (p instead enters via the the matter 
Lagrangian, which in the Einstein frame couples to both 
the metric and the scalar field, and shows up as a scaling 
of masses in the Einstein frame compared to masses in the 
Jordan frame. This scaling can be deduced from the rest 



mass of a particle: m 



F(ip)m 



2 

Jordan 



Einstein ' 



That is, 



^Einstein 



^Jordan 



F(ip)g^p^ 



(35) 



Hence, setting Mb = 1 corresponds to using TeV units in 
the Jordan frame, but Planck units in the Einstein frame. 

Using Eqn. I|35|l together with the usual scale-factor 
dependence of a 3-volume, one finds 




e -V+¥>o 



(36) 
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where R denotes the Einstein frame scale factor and Rq, 
ifo and A refer to the present day values of the quan- 
tities. In contrast to pressureless matter, the radiation 
density p r evolves normally (p r oc i?~ 4 ). This can be 
viewed to arise as a result of the zero rest mass of the 
photon. Including the contributions from x on e eventu- 
ally finds: 

P = fl-O^j-J +P <°{-r) 

+ \x 2 + V( X ) , (37) 

P = Pro^ 4 + lx 2 -V(x). (38) 

Equations (|20II22|) and (|26H28I) together with the Fried- 
mann equations I|29II30[) and the explicit expressions for 
density and pressure ((SHHSHIl, form a complete set of 
equations for the evolution of x an d the Einstein frame 
scale factor R. 



IV. STABILIZATION AND CONVERGENCE OF 
PERTURBATION THEORY 

Section Hi CI presented an argument for stabilization of 
a finite radion in the LED context, based on a particu- 
lar form of the Casimir potential, Eqn. 11411 . According 
to Eqn. l(T5|l the stabilization with e v = M h r ~ 10 15 is 
consistent with the perturbative scale |c| ~ e ~ 1/50. At 
first, such a small scale looks perfect for a perturbative 
treatment of the problem. This is good since the action 
JHJ is based on an effective one-loop parametrization of 
the radiative corrections. There is a fundamental flaw in 
the above argumentation however, since the actual ex- 
pansion parameter is not e, but 

e = eip . (39) 

Now, reading Eqn. I|15(l the other way around, we see 
that for tp = log Mb r » 1 it sets e ~ l/<f, or equivalently 
e ~ l! Hence, the perturbative expansion appears to be 
necessarily diverging at least in the case of the poten- 
tial. But since all corrections originate from the same 
bulk scalar field <\> in the ABRS scenario, they should 
all be roughly of same order a ^ b ^ c ^ 1/50, which 
renders the entire one-loop parametrization in Eqns. 
lllll useless. In fact the stabilization mechanism itself 
would only seem to work qualitatively at best, and could 
be destroyed by unknown higher order loop and non- 
pcrturbative corrections. 

It is quite interesting that a similar conclusion also fol- 
lows from the observational constraints on gravity. Pre- 
cision tests of General Relativity in the solar system sets 
a bound on the coupling between the scalar field x an d 
matter [l2|: 

a 2 ( X o) <2.2x 10- 4 = a§. ( 40 ) 



From Eqn. (|28|l it is then clear that one must have very 
accurately A ~ —a/2. From the one-loop expression A = 
1+ap and the fact that <p is large by assumption, we then 
see that a < and \a\ <C 1, so that a solution only exists 
if atp ~ —1. That is, the one-loop contribution must 
again balance the tree level one to a high precision. 

We started out by wanting to view the expansions in 
Eqns. l|9"l lll[) as coming from radiative corrections. How- 
ever, since the radiative corrections need to be large, they 
cast doubt over the self-consistency of the theory. While 
the precise forms of the expansions is ultimately irrele- 
vant - all that matters is that A has a zero at large ip 
(corresponding to a logarithmically large value of M\>r in 
the LED context) - it is interesting to see that Eqns. 
lll|l can be obtained in the context of a converging per- 
turbation theory. Quite simply, one only needs to as- 
sume that the radiative corrections arise not from one, 
but from several and possibly very many independent, 
or at most weakly coupled bulk fields (f>i. In such a case 
the loop corrections from each independent bulk field 4>i 
may be small, while their combined one-loop correction 
can add to a value close to unity. To illustrate the idea 
assume that there are N >■ 1 completely independent 
scalar fields <pi with equal cubic self-couplings living in 
the bulk. One would then have a truly small perturba- 
tion expansion parameter 

e N = e/N~l/N (41) 

for each separate bulk field In this case, the full per- 
turbation expansions for the coefficients A 7 B and C ap- 
pearing in the action (jSJ will be converging. Schemati- 
cally for example: 

N 

A = Y, A = 1 + #Ne N + #iV4 + . . . , (42) 

i=l 

where #:s represent numbers of order unity. If we now as- 
sume that the one-loop contribution is of order one, then 
the higher order loop corrections are suppressed by pow- 
ers of 1/iV; for example at the two-loop Ne 2 N ~ e/N ~ 
1/N. Of course it is not necessary to assume that all 
bulk fields <j>i have exactly the same couplings. It is 
enough that the sum of radiative corrections is chosen 
such that the combined one-loop contribution is of order 
one at the present value of if. On the other hand, one 
must prevent large mixing couplings gijk in order to keep 
higher order loop corrections small. Indeed, at two-loop 
in the strongly couplicd limit one would find for example 
5A& ~ N 2 e 2 N ~ 1. 

The same qualitative argument also apply to the func- 
tion B. Tuning A close to zero does not imply that B is 
small at the same time however. Still, the neglected two- 
loop corrections are small for both quantities. The case 
of the Casimir potential in the LED context is different; 
for Uc there is no tree level contribution as the whole po- 
tential is purely of quantum origin. Its multi-field version 
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corresponding to Eqn. (|42(l then looks as follows: 

N 

Uc = ^U i = NU N {l + #e N + #e 2 N + ...), (43) 

i=l 

After we make the obvious identification NUn — Uq ~ 1 
we observe that the radiative correctios to Uq are truly 
small: c ~ #£jv ~ 1/N. This implies that the original 
ABRS-type stabilization relying on the corrected Casimir 
potential is not feasible in this construction. 

A. The New Stabilization Mechanism 

Fortunately, the model does in fact have another sta- 
bilization mechanism built in it. Indeed, if the one-loop 
correction, #Nsn = atp in Eqn. (|42() is negative, then the 
function A(tp) will become zero in the scale tp ~ 1/ |a| . 
Since the two-loop corrections to A, B and C are still 
small at this point, they can not alter this conclusion 
significantly. The new stabilization mechanism now sim- 
ply follows from the l/^4 2 -scaling of Ve(<^), which makes 
it an increasing function of tp as A approaches zero, and 
eventually gives rise to a narrow but infinitely high con- 
fining potential barrier at A = 0. 

The potential shown in Fig. Qd) actually corresponds 
to A — confinement. The minimum seen there has 
nothing to do with the potential U , for which we have in 
fact used a pure exponential, runaway form with c = 0. 
Moreover, the sharp feature at the top of the poten- 
tial maximum, the kink, corresponds to the singularity 
A(x*) = 0) which would have shown up as an infintcly 
high and narrow spike in Ve (</?). The fact that the singu- 
larity becomes a finite height kink in the x-picture follows 
from the highly nontrivial relation between \ and tp (see 
appendix A). One can understand this when one observes 
that also k 2 {tp), Eqn. 1)20(1 . becomes singular at A = 0, 
so that also the kinetic term for tp diverges. This has the 
opposite sign in the action and balances out the potential 
singularity. (Obviously the model should be qualified as a 
mixed fc-essence and quintessence in the singular region.) 
Although V(x*) is finite, note that it still is 20 orders of 
magnitude larger than the value of the potential at the 
minimum! Together with Hubble friction this barrier be- 
comes impenetrable, prohibiting the field from reaching 
values beyond \* an d thus giving rise to a robust con- 
finement mechanism for tp, corresponding to large, finite 
radii in the LED context. 

In summary, we have shown that in a scenario with 
many weakly coupled independent bulk fields the ef- 
fective theory Eqn. © can be reliable even when A, B ss 
0. Stabilization for large values of tp. corresponding to 
large radius in the LED context, can not be obtained 
from the corrected Casimir potential because the scenario 
predicts C ~ 1. However, it can originate from a singu- 
larity in the potential at A = 0. On the technical side, 
the fact that the full one-loop contributions to quantities 
A and B are large, implies that one can not simplify any 



of the expressions that depend on them in the previous 
sections, by expanding in cap or bip as was done in ref. 
This is of course the reason why we always kept the exact 
expressions. Finally, it should be stressed that while our 
quantitative results will depend on the particular form 
of the one-loop parametrization, the qualitative features 
would be the same for any function A that crosses zero 
at a sufficiently large values of tp. 



V. CONSTRAINTS 

The underlying model JSJl has four unknown param- 
eters: £/o, a, b, and c. Given a full knowledge of the 
physics leading to the effective theory, these parameters 
could be calculated exactly in terms of the fundamental 
coupling constants, but as those details are not known, 
one has to settle with combinations of phenomenological 
and theoretical naturalness constraints. Stabilization ar- 
guments, together with the desired large values of r in 
the LED context, have already settled the overall scale 
of the parameters Uq, a, b and c. We shall next impose 
the constraints coming from various observations. 



A. Observational Constraints 

There are four major observational constraints affect- 
ing the model parameters. We have already mentioned 
the gravity constraint on the scalar field coupling a at 
present, and we shall return to it again shortly. The 
other constraints are of astrophysical origin, mainly orig- 
inating in the CMBR and supernovae observations. As 
previously stated, it is in the Jordan frame where all 
observations have their usual interpretation. Hence, to 
obtain the true constraints on the model coming from for 
instance the CMBR, one needs to make a full computa- 
tion of the possible spectra coming from the model and 
fit them to the observed data. This is however beyond 
the scope of the present paper. 

Instead, we will use a more modest approach based 
on the fact that the behaviour of our model eventually 
is very close to that of General Relativity. In such a 
case, F{tp) and Z(tp) are essentially constant so that Jor- 
dan and Einstein frame quantities do not differ except for 
overall constant factors. E.g., given that M 2 F(tp) ~ Mpj, 
setting Mb to one corresponds to using TeV units in the 
Jordan frame, but Planck units in the Einstein frame. 
The majority of the constraints concerns the present day 
value of various quantities. Hence, applying the con- 
straints will essentially only give us the value of our pa- 
rameters at one specific point (the present) in the evolu- 
tion of the universe. In practice, we then start the evolu- 
tion of the universe and shoot for this point. Given that 
we succeed, we then need to check that the behaviour of 
the model was close enough to General Realtivity at ear- 
lier times. Let us now return to discuss these contraints 
case by case. 
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The dark energy component is constrained by it's 
present amount and by the equation-of-state parame- 
ter w x . The equation-of-state parameter is limited by 
CMBR and astronomical data to 

w Xo < -0.78 . (44) 

To be precise, the limit (|44f) applies for the present value 
of w x , and was derived in standard General Relativ- 
ity with the assumption that both w x and p x are con- 
stants. In the present model both the energy density 
and the pressure of the quintessence field x are dynam- 
ical variables and so the bound is not directly ap- 
plicable. Moreover, the field x a ^ so couples to regular 
matter. More general limits on the function w x (z) have 
been found in the literature j2^,|23. However, we take 
Eqn. I|44|l as an indication that the present value of w x 
should be very close to —1. It is convenient to rephrase 
Eqn. 144(1 in terms of present ratio of the kinetic and the 
potential energy densities of the field x : 

S X0 = \^ < 0-12 ■ (45) 

The third observational constraint comes from con- 
necting the model parameters with the observed amount 
of dark energy at present: 

V( X o) = MtU %-^ = J^L, (46) 

A 1 + X0 

where J7a ~ 0.73. The content of Eqn. (|46|l is that a given 
set of model parameters Uq, a and c, and the boundary 
condition for S Xo will determine the precise value of tpo, 
i.e the value of x today. Because p c ~ 10~ 120 M^ in the 
Einstein frame, it is clear that in the LED context, setting 
Uo ~ 1 imposes roughly the desired scale of M^r — e v ~ 
10 15 . 

Finally, there is a constraint due to the need of preserv- 
ing the success of big bang nucleosynthesis: The universe 
must be radiation dominated already at BBN-time, i.e. 
when the observed Jordan frame temperature Tj ~ 1 
MeV. This constrains both the kinetic and potential en- 
ergy densities of the scalar field, which would otherwise 
make the universe expand too fast and result in too large 
helium production. This consideration does not lead to 
additional constraints on model parameters, but rather 
to a need for tuning the initial conditions for the field. 
There is a correlation between the two however, because 
different sets of model parameters need different sets of 
initial conditions to provide an acceptable cosmological 
solution. Moreover, it is not a priori clear that accept- 
able solutions do exist at all for arbitrary sets of model 
parameters. Obviously these questions can ultimately be 
studied only through numerical analysis. 

B. Fine- Tuning in the Model Parameter Space 

Let us now return to the gravity constraint Eqn. I|40|) 
(a-bound for short). While the form of the potential 
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FIG. 2: Shown is the present value of the square of the cou- 
pling a as a function of the parameter b for three different 
values of the parameter a. The other model parameters are 
set to the values given in Eqn. 1471 . The a-bound, Eqn. 1401 . 
allows the region below the dotted line. 

automatically adjusts A close to zero, it is not immiedi- 
ately clear whether this tuning is sufficient to satisfy the 
a-bound automatically. Comparing Eqns. i|27|) and l|28|) 
we see that there is in fact some tension between fulfilling 
the a-bound and reaching the potential minimum where 
V(x) = 0- Indeed, for c = the latter corresponds to 
A = —a/4, whereas the former goes to zero at A = —a/2. 
Whether this difference can be accommodated obviously 
depends on how strict the a-bound really is. From the 
parametric dependence of a = a(x, a, 6) one sees that the 
a-bound leads to correlations between the model param- 
eters, and so does the need to adjust V — V(x, Uo, a, c) 
to the present cosmological constant as well. In the end 
one is forced to some cross-correlated fine-tuning with 
the model parameters. Instead of trying to map out the 
allowed region in the parameter space, we estimate the 
degree of the fine-tuning that is needed by fixing all other 
parameters to some natural values and then checking how 
much freedom is left in varying a. Observe that this does 
not mean that a is a priori the parameter to be fine 
tuned; we could as well set a rather freely, and then fine- 
tune either Uq, or even xo, i-e. the radius of the extra 
dimensions today ro in the LED context. 

Let us assume that xo is fixed. For definiteness we 
assume that the potential parameters and 6 Xo take on 
the values 

Uo = 2.60 , c = and (5 Xo = 10~ 2 . (47) 

For a fixed xo the a-bound implies that Aq + a/2 has to 
be fine-tuned to zero to an accuracy 

\A + a/2\ < ^Ja 2 \aB \/2. (48) 

Because of the smallness of a§ this is always a rather 
strong constraint, and it becomes particularily severe 
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FIG. 3: k 2 (ip) as a function of ip/ Tog 10 (= log 10 M^r) for different values of the parameters a and b. Thick gray lines highlight 
the values of ip corresponding to a phantom phase. Shaded regions correspond to values of ip where gravitons would carry 
negative energy in the Jordan frame. 



when | So | is small. This dependence is clearly visible 
in Fig. J2J), where we plot a as a function of b for some 
representative values of a. For the values chosen above, 
l-Bol goes to zero at « —0.0280 0] and the closer b 
gets to &*, the more accurately a must be fine tuned to 
a value a — » 6*/(l — &*/2) —0.0276. If we now impose 
the perturbative argument that \b\ is of the same order as 
\a\, it implies that a has to be fine-tuned to roughly one 
part in ten to hundred thousand, which corresponds to 
a relative accuracy of somewhat less than one per cent. 
This gives an estimate of the amount of fine-tuning that 
is needed when adjusting different acceptable model pa- 
rameter sets. In what follows we shall fix a to the value 
that allows the widest interval in b in Fig. J5J), i.e. 

a = -0.02759 and b < -0.0300 . (49) 

Eqns. H47f) and (|49() define the allowed parameter space 
for our quintessence-like scenarios with the model ac- 
tion Changing Uq and c only moves solutions to 
a slightly different band in a and b space. To be specific, 
we will choose b = —0.034 as our reference value in what 
follows. 

Finally note that the expression on the right-hand side 
of Eqn. lf4*%|l in general is smaller than |a|/4 » 0.007, the 
value of |A| corresponding to the minimum of V(x)- I n ~ 
deed, for our reference value of b we have y a§|a_Bo| ~ 
0.001. As a result the cosmological evolution must be 
such that the field is x actually not in its potential mini- 
mum at present. Moreover, all our solutions should show 
some late time evolution and adjustment of a close to 
the present time. One can lessen this effect by choos- 
ing a large value for \Bq\, but apart from such additional 
tuning the late evolution should be a generic prediction 
of the model. 



C. Phantom Parametrizations 

For our reference set, Eqns. (|47() and l|49|) with b = 
—0.034, the model stays non-phantom for all values of ip 
with F(ip) > 0. As seen in Fig. (PK), the field ip becomes 
phantom over a finite interval in the region where gravi- 
tons would carry negative energy in the Jordan frame, 
but since the field evolution never enters the this region, 
such cosmologies cannot be realized from natural initial 
conditions. However, it is also possible to arrange the 
model parameters such that <p becomes a phantom field 
already before the A = singularity. We show examples 
of such models in Fig. (J3Jd). It is quite interesting that 
these phantom-like structures were found using perfectly 
natural values for the model parameters. Nevertheless, 
it is not clear to us whether one can arrange a phan- 
tom phase into the past, present or the future of our 
universe such that it would be compatible with all the 
observational constraints. At least the examples that we 
have studied connected to large extra dimensions, includ- 
ing the models shown in Fig. J3t>), almost certainly are 
ruled out by gravity constraints. However, if one indeed 
takes the generic scalar-tensor theory of Eqn. as the 
starting point, several restrictions such as the millime- 
ter scale for large extra dimensions no longer apply, and 
one could obtain a richer structure of possibly acceptable 
mixed phantom-quintessence-fc-essence solutions. 

A promising study addressing these questions has been 
done by Perivolaropoulos [3l| who show that, in contrast 
to the case of a minimally coupled scalar field [32], it is 
indeed possible to reproduce an eventual phantom transi- 
tion within scalar-tensor theories. Perivolaropoulos em- 
ploys a link between the equation of state of dark energy 
and the precision tests of General Relativity in the so- 
lar system. In this scheme, one redefines ip in the Jordan 
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TABLE I: Initial conditions/starting values and the resulting 
values at present, given in Planck units. Note that p to t = p c - 

frame so that the field has a canonical kinetic term (apart 
from the sign). Supernovae data and the GR-constraint 
are then used to put bounds on the redshift z dependence 
of (f, the potential U, and the function F multiplying R g 
in the Jordan frame action. By assuming a particular 
form of ip{z) it is then possible to reconstruct U(ip) and 
F(ip). We leave such explorations for upcoming studies. 
The conditions for obtaining a dark energy equation of 
state less than —1 has also been explored in ref. [3^|. 



D. Other Constraints 

Some of the constraints discussed so far were generic 
to any effective quintessence model. Others, such as the 
need for ip to have the same value at BBN-scale as at 
present, were particular to the present model having its 
origin in scalar-tensor theory. In the LED context, there 
is yet another type of constraint particular to these mod- 
els that we need to discuss. Namely, the radion degree 
of freedom that we have seen to give rise to quintessence 
field, comes associated with an infinite tower of Kaluza- 
Klein excitations on the brane. The problem with the 
KK-modes is that if the universe enters the radiation 
dominance at too high temperature T*, too many KK 
modes might be excited through their couplings with or- 
dinary particles [2(], 0, H^- However, it is quite diffi- 
cult to state precisely how these considerations limit the 
model parameters, because they depend on introducing 
an entire new structure to the model, which has not yet 
been explored in detail. Indeed, the problem of KK- 
modes might be associated and solved by the unknown 
higher-dimensional physics, in which case even no ad- 
ditional constraints might arise @- Even in the more 
conservative view, which assumes the possiblity of ex- 
citing the KK-modes at the reheating scale T*, bounds 
depend on the details of the temperature history on the 
brane. They can be quite stringent however; for example, 
Hannestad pjH ] reports that if the universe enters radia- 
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FIG. 4: The relative difference between the Einstein and the 
Jordan frame scale factors AR = (R—Rj)/R (solid) and times 
At = (t — tj)/t (dashed), where we have set R = Rj = 1 and 
t — tj at present. 

tion dominace instantly, then the fundamental mass scale 
(having fixed M h r) is bounded by M h > 8 + lO^/MeV) 
TeV. In principle we can always accommodate such a con- 
straint in our model. This would of course to some degree 
undermine the particle physics motivated attraction to it, 
which originates from connecting the fundamental grav- 
ity and electroweak scales to evade the heirarchy prob- 
lem. From a cosmology point of view, the main problem 
is however that an increase in Mb must be compensated 
by a decrease in r, which would eventually spoil the natu- 
ralness of the Casimir potential as the source of dark en- 
ergy. The situation is not extreme however; taking a re- 
heating temperature T„ = 1 MeV, gives r < (0.02eV) _1 , 
which would still imply a reasonably modest fine-tuning 
of order Uq ~ 10~ 5 . 

VI. A REALISTIC COSMOLOGY 

In the previous section we showed that setting the fun- 
damental energy scale to TeV, or equivalently, the stabi- 
lized extra dimensions to millimeter scale in the LED con- 
text, combined with theoretical and observational con- 
straints practically fixes the fundamental parameters in 
the model. This means that the form of the quintessence 
potential V(x) is almost set as well, not leaving much ex- 
tra space for adjustments to get an acceptable cosmolog- 
ical evolution. Of course, a lot has already been achieved 
in the course of fixing the parameters; built in is the fea- 
ture that the potential has a minimum, corresponding 
to millimeter scale in the LED context, such that V(x) 
would give the right amount of dark energy today. Also, 
the gravity constraints and constraints on the equation- 
of-state parameter are almost automatically satisfied if 
the field sits near the potential minimum. However, we 
still have to show that the model indeed does behave 
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FIG. 5: Different components of the energy density (in Planck 
units) as a function of the Einstein frame scale factor R. 
The various curves represent matter (dashed black), radia- 
tion (dotted) and scalar energy density (solid). The dashed 
grey and dot-dashed curves shows the potential and the ki- 
netic part of the scalar energy density, respectively. 



very similar to General Relativity and that it can support 
a cosmological evolution consistent with all constraints, 
leading to the desired densities at present. 

Let us stress that we actually do not know the full 
theory at the fundamental TcV-scalc; we have to our dis- 
posal only an effective theory that is valid at much lower 
energy scales. In fact, the thermal history beyond the 
entrance to radiation dominance is subject to unknown 
physics such as Kaluza-Klein excitations in the LED con- 
text and the constraints arising thereof PHHH. We 
will therefore start off the evolution of our model universe 
at a temperature well below the TeV-scale, Tj <§; 1 TeV, 
but still above the big bang nucleosynthesis temperature 
of Tj ~ 1 MeV. In practice we are using the Einstein 
frame scale factor R as an integration variable and we set 
R now = Ri ow ee 1. We know that R^ BN ~ 10" 10 - 1(T 9 , 
and i?xeV — 10~ 16 , but it is not possible to determine 
the corresponding values of the Einstein frame scale fac- 
tor until we know the actual evolution. However, since we 
have assumed that the model should behave very close to 
General Relativity, the Einstein frame scale factor should 
not differ significantly from the Jordan frame scale fac- 
tor. Having chosen a suitable starting point i? s tart, we 
integrate the evolution equations until present with vari- 
able initial conditions for \ and \ until (if) an acceptable 
solution is found. In this scheme, the first and most im- 
portant test is to see if our evolution indeed behaves like 
General Relativity after nucleosynthesis. Having estab- 
lished that, we can then check if other quantities such as 
energy densities fulfils the desired constraints. 

We now turn to the qualitative features to be met by 
an acceptable solution. First one has to make sure that 
the model does not spoil the BBN explanation for the 
observed light element abundances. As was explained in 



FIG. 6: The equation-of-state parameters w x (solid) and uitot 
(dashed) as functions of the Einstein frame scale factor R. 



the previous section, this is ensured if tp was already very 
close to its present value during BBN. This turns out to 
be a tough constraint, which basically forces us to set the 
field relatively close to its potential minimum at i? s tart- 
Fixing BBN combined with the constraints built into our 
model potential, go a long way towards fixing the correct 
cosmological model today. What remains is essentially 
the issue of getting the quantitative values of present en- 
ergy densities J7 X ~ 0.73 and Sl m w 0.27 to coincide with 
Ptot ~ Pc- This of course is just the familiar "why now" 
problem, which in general requires some fine-tuning of 
initial conditions in generic quintessence models. Here 
the "why now" problem is essentially transformed to the 
assumed TeV scale of Mb, corresponding to millimeter 
size extra dimensions in the LED context. 

Below we present an example of a realistic cosmol- 
ogy corresponding to our basis set of model parameters; 
U = 2.60, a = -0.02759, b = -0.034, and c = 0. 
We have started the evolution at i? s tart = 10~ 12 , which 
roughly corresponds to the scale T s J tart ~ 100 MeV. The 
corresponding initial conditions/present values of x, x 
etc. are given in Table HJ and the cosmological evolution 
of Rj and ij, energy densities pi, the equation-of-state 
parameter w x , the field tp and the coupling a are shown 
in Figs. 10} , 101, ©, OO) and Ct>), respectively. Note 
that we have extended the evolution a bit further than 
present, with grey areas representing the future. 

First of all, we see from Fig. that, except at very 
early times, the relative differences between the Einstein 
and the Jordan frame times and scale factors are neglible, 
so that the frames are descriptively equivalent. That is, 
the rate of Einstein frame clocks and size of Einstein 
frame rods essentially do not differ from the standard 
Jordan frame clocks and rods except for a constant scal- 
ing. Hence, our model does in practice behave like Gen- 
eral Relativity and our application of the constraints in 
section IV Al remain valid. This also tells us that ip has 
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FIG. 7: (a) Shown is p/log 10 (= log 10 M^r) (solid line) and the fc 2 (ip)-factor defined in Eqn. 11201 (dotted line) as a function 
of the Einstein frame scale factor R. (b) The coupling a as a function of the Einstein frame scale factor R. 



remained very close to constant, so that the strength of 
gravity had the same value at BBN as at present. Let us 
now proceed with discussing the evolution in the Einstein 
frame in detail. 

Due to the relatively large starting value of ip, the po- 
tential is neglible compared to the kinetic energy and the 
scalar field x enters radiation domination in a kinetic- 
dominated roll. During early stages, the kinetic energy 
scales as \ 2 ^ R~ 6 , and it dominates the energy density 
for log 10 R < —11. This scaling follows from Eqn. i|26|) 
when the potential and the a-term are neglected and the 
evolution of the scalar field is dominated by the Hubble 
friction 3Hx- During this kination phase |3(| the field 
X is rolling down its potential and ip, corresponding to 
the radius in the LED context, is growing rapidly. The 
scalar field dominance ends and the universe enters radi- 
ation domination right before BBN. At this time x settles 
close to the minimum of the potential, and correspond- 
ingly ip reaches a plateau as shown in Fig. $7\i). The 
energy of the scalar field continues to fall rapidly, un- 
til kination is halted by the a-term. The a-term causes 
X 2 to track matter density until log 10 R « —1.5, corre- 
sponding to and observed redshift zj « 30 0], where 
the kinetic energy falls below the potential, and the field 
X starts acting like a cosmological constant. 

Although V{x) now dominates the total energy den- 
sity, it is not yet a constant (see the linear scale plot in- 
serted in Fig. 10). The cause for this final 10% decrease 
between zj « 0.5 and present can be tracked to the ten- 
sion between the a-bound and the position of the poten- 
tial minimum for our model parameters, discussed in the 
last paragraph of section IV Bl Indeed, for a long time, 
the a-term acts like an additional effective potential in 
the equation of motion l|26(l and confines the field to the 
proximity of the zero of the a- function at A ~ —a/2. 
As the matter density becomes very small, the a-term 
eventually becomes small in comparision to the poten- 



tial gradient, and the field can finally roll down to its 
true minimum. This phase can also be seen as a slight 
increase in ip in Fig. GK), as well as in a little more dra- 
matic change in k 2 (tp) and finally the coupling a itself 
shown in Fig. Qd). 

The evolution of both w x and w to t are shown in 
Fig. 0. The epochs of kination (w to t = 1), radiation 
dominance (wtot = 1/3), matter dominance (wtot = 0), 
and the eventual cosmological constant dominance when 
Wtot — ► — 1 are clearly visible. The evolution of w x is 
much less dramatic. It remains close to unity all the way 
until log 10 R ss — 2, when it rapidly falls from w x = 1 to 
close to w x = — 1, as the kinetic energy of the field be- 
comes neglible in comparision to the potential energy. 
The only interesting feature is the little bump in w x 
which is caused by the small temporary increase in the 
kinetic energy of p x , as % makes its last roll closer to the 
potential minimum. 

In summary, it is clear from Figs. I@J and and Ta- 
ble |U that, except at very early times, the Einstein and 
Jordan frame are descriptively equivalent and the evo- 
lution of all densities have the desired features. In par- 
ticular, the universe is flat with f2 x = 0.73 at present. 
Also, as shown in Figs. (JTJt) and (0)), both ip and a 
fulfill the observational constraints. Overall, the evolu- 
tion of the energy densities is qualitatively surprisingly 
similar to the results of Albrecht et al. 8), apart from 
the additional evolution at late times connected with the 
mismatch between the zero of the coupling a(x) l|28|) and 
the minimum of the potential V(x)- 



A. Quintessence Vs. A-Cosmology 

Although the scalar field behaves like a cosmological 
constant at present, it obviously had a varying equa- 
tion of state in the past. This allows for a possibility to 
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FIG. 8: (a) Observed redshift-magnitude diagram normalized to a flat, J1a = 0.73 cosmology, with binned simulated SNAP 
data fl^| . The solid line representing the quintessence-like scenario corresponds to Sl x = 0.73 at present, (b) Statistical X 2 (°"i) 
for a cosmological constant, using the quintessence-like scenario with Q x = 0.73 at present as null hypothesis. 



put bounds on the model based on observations that can 
probe the equation of state at different times. Moreover, 
the proper way to confront a scalar-tensor theory with 
supernovae observations is to compute its corresponding 
magnitude-redshift relation. This focus on the observed 
expansion rate alone and does not depend on some par- 
ticular parametrization of the theory. We have therefore 
computed the magnitude-redshift relation of type la su- 
pernovae, predicted by our reference case. The result is 
given in Fig. (|BK)j where the magnitude-redshift relation 
is shown normalized to a corresponding flat, £l\ = 0.73 
cosmology. The model easily fits current data £3] ■ How- 
ever, the magnitude-redshift relation will be measured to 
high accuracy by the SNAP-satellite 13] in the near fu- 
ture. Comparing with the shown binned simulated SNAP 
data obtained from ref . [l4| , suggests that the satellite ex- 
periment will be able to distinguish the models from each 
other. This interpretation is verified by Fig. JHJd), which 
shows a X 2 (<7i) function in the space of constant A mod- 
els using our quintessence-like scenario with Q Xo = 0.73 
as null hypothesis. The f^A = 0.73 model lies clearly 
several sigmas away from the present model. However, 
the minimum of the % 2 -function at Q\ = 0.70 lies well 
within one sigma from the null hypothesis. Thus, SNAP 
data alone will not able to distinguish between our model 
and the A models with about 5% less dark energy. While 
the difference is small, it might be possible to break the 
degeneracy by combining SNAP-data with current data 
on the cosmic microwave background anisotropies. 

Let us finally remember that the late time evolution, 
due to the tension between the zero of the coupling a and 
the minimum of the potential V, is a rather generic fea- 
ture of the present model, and therefore so is the observed 
difference in the magnitude-redshift relation between the 
present model and a corresponding constant A model. 



B. Robustness of the New Stabilization Mechanism 

Finally, let us consider the robustness of the new sta- 
bilization mechanism based on the A = singularity. On 
one hand, the barrier height at the singularity is tiny 
compared to the value of the potential when tp w Q, but 
on the other, it is about 20 orders of magnitude higher 
than the value of the potential at the minimum. Could it 
nevertheless be possible to go past the barrier and make 
the world enter the region where gravitons would carry 
negative energy in the Jordan frame by giving the field 
a much larger initial velocity? The answer turns out to 
be negative, for all initial velocity we gave the field, the 
solution always either got stuck close to the singularity, 
or got its kinetic energy diluted by the Hubble friction. 

Figs. (JSJi m>) shows an example of an evolution with 
huge initial field velocity, where we have used the same 
reference values of the model parameters and initial con- 
ditions as shown in Table [fl apart from giving the field 
12 orders of magnitude larger kinetic energy. Now the 
potential reaches its minimum already at log 10 R « — 1 1 
and starts increasing again as the field continues climb- 
ing the potential barrier. The growth of x an d therefore 
also that of tp stops at log 10 R w — 6, where the potential 
energy reaches its maximum and the kinetic energy goes 
through zero (this shows up as a small dip in Fig. JHK))- 
The huge spike in the A; 2 ((^-function shows that the field 
came close to the singularity before turning back. As the 
field rolls back, ip is changing slowly initially, because of 
the near singular relation between % and tp, but even- 
tually at about log 10 R ~ 3 (far in the future) tp makes 
a small adjustment before freezing out completely. The 
field then rolls back towards the minimum as the Hubble 
friction absorbs the remaining kinetic energy. Observe 
that most of the time during the period of backwards 
rolling the universe is neither matter, radiation, or cos- 



16 



Fig. (a) 



Fig. (b) 



-60 



-80 



° -100 
o 



-120 



-140 



-160 



\ ' • , N. 
















s 

















15.5 



-12 



-2 3 
logio R 




14.5 



FIG. 9: (a) Shown are energy densities (in Planck units) as a function of the Einstein frame scale factor R for the same model 
parameters and initial conditions as in Fig. ||KJ, except for the initial field velocity which now has been set to dx/dt = 2.0 x 10 -31 
Ml. Different curves represent matter (dashed black), radiation (dotted), total scalar energy density (solid), scalar potential 
energy (grey dashed) and scalar kinetic energy (dot-dashed), (b) Shown is ip/ log 10 (= log 10 M^r) (solid line) and the k 2 (ip)- 
factor defined in Eqn. 1201 (dotted line) as a function of the Einstein frame scale factor R. 



mological constant dominated, nor does it behave like 
during x 2 -dominated kination. Instead, it is a peculiar 
mix of kinetic (fc-essence) and potential (quintessence) 
energy. Needless to say, this example is entirely incom- 
patible with our universe; we only gave it to show the 
robustness of the radius stabilization against initial con- 
ditions with very large field velocities. 



VII. SUMMARY AND DISCUSSION 

We have considered in detail the cosmologies of a class 
of dilatonic scalar-tensor theories of gravity inspired by 
large extra dimensions and based on ideas put forward 
by Albrecht, Burgess, Ravndall and Skordis 0, 0j- ^ n 
the LED context, the effective quintessence field \ cor ~ 
responds to a radion mode originating from dimensional 
reduction of a full (4 + d)-dimensional action to an effec- 
tive 4D brane action. We have shown that observational 
gravity constraints (the coupling a) require that radiative 
corrections from additional bulk fields 4»i nearly cancel 
the tree level contribution to the effective gravitational 
constant G e g at present. The corresponding zero in the 
scaling function A(ip) of the Ricci scalar in the Jordan 
frame, shows up as a singularity in the effective potential 
in the Einstein frame, yielding a new type of stabilization 
for the field ip (i.e. the radius in the LED context). In 
this stabilization scheme the model is thus (almost) auto- 
matically consistent with the gravity constraints. If this 
mechanism is implemented in a perturbatively converg- 
ing scheme, the lowest order corrections to the Casimir 
potential are truly small, and can not give rise to a sta- 
bilization, as was assumed in the original work 0|. 

Depending on the parameters, the model can give rise 
to features very different from a generic quintessence 



model. To a degree this was already recognized by the 
earlier work, and indeed the model has been dubbed 
"walking quintessence" 0- However, we have shown that 
the model can behave at times more like /c-essence, where 
the dynamics is mainly controlled by the kinetic term, or 
even as a phantom cosmology, where the scalar field has 
negative kinetic energy. Using perfectly natural values 
of the model parameters, it is even possible to arrange a 
phantom phase that has finite range in ip, which might 
allow for a finite duration phantom phase, possibly side- 
stepping the associated problems with vacuum instability 
in such models However, whether such a region can 
be reached dynamically in a model compatible with the 
observational constraints is yet unknown to us; although 
it seems unlikely in the subset of models studied in this 
paper, the possibility at large remains an open question. 

Imposing various constraints very nearly determine all 
the fundamental parameters in the effective 4D brane ac- 
tion. Moreover, some of the parameters or combinations 
of parameters needs to be fine tuned for the model to 
work. (That is, to simultaneously give rise to the ob- 
served strength of gravity and obey the a-bound.) Nev- 
ertheless, having fixed the parameters, we show that the 
model can give rise to an acceptable cosmological evo- 
lution. Moreover, the model predicts a characteristic 
signature in the magnitude-redshift relation of type la 
supernovae, that may make it possible for a combination 
of supernova and CMBR data to distinguish the model 
from a pure A-cosmology. When applying the observa- 
tional constraints, we made the assumption that the be- 
haviour of our model should be very close to General 
Relativity. Although beyond the scope of this paper, a 
complete observational analysis requires a computation 
of the corresponding CMBR spectrum and the cosmolog- 
ical perturbations constrained by large scale structures. 
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A particularily attractive feature of the present model 
is that it makes a clear identification for the origin of the 
effective quintessence field in a well defined context that 
has other motivations beyond the cosmological applica- 
tion. Indeed, models with large extra dimensions were 
originall y propo sed as a novel solution to the heirarchy 
problem |20ll2ll |. This blessing comes with a burden how- 
ever. We have seen that cosmology puts several stringent 
constraints on the model parameters and leaves little free- 
dom in choosing the initial conditions of the cosmologi- 
cal evolution. A particularily unappealing feature in this 
context is that one is forced to assume a very large value 
for the radius r at the time when the universe enters 
the radiation dominance. Indeed, the thermal history in 
LED models is heavily constrained by the need to avoid 
overproducing Kaluza-Klein modes, and not much can be 
said with a certainty about the conditions prevailing at 
or before the reheating phase. We did not discuss these 
issues in any detail in this work, and refer the reader to 
the original papers 0, @ for elaboration in this matter. 

While it is certainly rewarding to be able to put the 
model into a larger theoretical framework, the scalar- 
tensor theory of Eqn. @ can indeed be viewed as the 
starting point. In this context the parametrization JSJ- 
lllll emerges as a simple correction to the usual dila- 
tonic scalar-tensor gravity. Adopting this point of view 
the model becomes much less severely constrained, and 
one could try also more complicated functional forms for 
the scaling functions A, B and C, for example to see if 
one could more easily construct phantom-like cosmolo- 
gies consistent with all astrophysical observations. 
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If j3 = the solution is instead 
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for 7 = one finds 



APPENDIX A: EXPLICIT RELATION 
BETWEEN x AND ip 



The defining relation for the field x is: 
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where A = 1 + cup, B = 1 + b<p and C = 1 + cip. The 
function appearing in the square brackets can be written 
as a polynominal in A = A/a: 
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x sgn A , (A9) 



and finally for A = 0, 
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All of the above solutions are valid in the range —0.1 < 
a < 0.1 and —0.1 < b < 0.1, given that they aand &are 
arranged so that k 2 {ip) > 0. 
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